Thermalisation time and specific heat of neutron stars crust 
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We discuss the thermalisation process of the neutron stars crust described by solving the heat 
transport equation with a microscopic input for the specific heat of baryonic matter. The heat 
equation is solved with initial conditions specific to a rapid cooling of the core. To calculate the 
specific heat of inner crust baryonic matter, i.e., nuclear clusters and unbound neutrons, we use the 
quasiparticle spectrum provided by the Hartree-Fock-Bogoliubov approach at finite temperature. 
In this framework we analyse the dependence of the crust thermalisation on pairing properties and 
on cluster structure of inner crust matter. It is shown that the pairing correlations reduce the crust 
thermalisation time by a very large fraction. The calculations show also that the nuclear clusters 
have a non-negligible influence on the time evolution of the surface temperature of the neutron star. 



I. INTRODUCTION 



The thermalisation process of the neutron stars crust 
can give important informations about the properties of 
the crust matter. This is indeed the case in the rapid 
cooling of isolated neutron stars [l[ and in the thermal af- 
terburst relaxation of neutron stars from X-ray transients 
In the rapid cooling process an important quantity 
is the cooling or thermalisation time of the crust, defined 
as the time needed for the crust matter to arrive to the 
temperature of the cold core. Since the outer crust and 
the envelope have high thermal conductivity, the cooling 
time is essentially determined by the inner crust mat- 
ter, formed by nuclear clusters, unbound neutrons and 
ultrarelativistic electrons. 

The thermalisation time of the crust depends essen- 
tially on the crust thickness However, several studies 
have shown that the thermalisation time depends also 
significantly on the superfluid properties of the inner 
crust baryonic matter [E y, 0, Q. This dependence is 
induced through the specific heat of unbound neutrons, 
strongly affected by the pairing energy gap. Since the 
neutron pairing gap is influenced by the presence of the 
nuclear clusters [3|, M 0| , a reliable calculation of neutron 
specific heat should take the clusters into account. How 
the intensity of pairing correlations affects the specific 
heat of the neutrons in the presence of nuclear clusters 
was analysed in Ref. Q. Thus, using the framework of 
Hartree-Fock-Bogoliubov (HFB) approach at finite tem- 
perature, it was shown that the specific heat can change 
with several orders of magnitude if the pairing gap is ad- 
justed to describe two possible scenarios for neutron mat- 
ter superfluidity, i.e., one corresponding to BCS approxi- 
mation and the other to calculations schemes which take 
into account in- medium effects [9(. The impact which 
these changes in the specific heat could have on the ther- 



malisation time was discussed in Ref. (4[. Em ploy ing 
a simple random walk model for the cooling (3. Ilfjj]. in 
which the diffusion of the heat towards the core was 
calculated without taking into account the dynamical 
change of the temperature through the whole crust, it 
was shown that the thermalisation times corresponding 
to the two pairing scenarios mentioned above differs by 
a large fraction, which could be easily discriminated ob- 
servationally. The scope of this paper is to present more 
accurate estimations of the thermalisation time obtained 
by employing a more realistic cooling model based on dy- 
namical solutions of the heat equations and on a state-of- 
the-art description of the specific heat for baryonic mat- 
ter. 



II. THE MODEL OF CRUST 
THERMALISATION 

The crust thermalisation is described here in the rapid 
cooling scenario in which the core arrives quickly to a 
much smaller temperature than the crust. Due to this 
temperature inversion the heat stored in the crust diffuses 
into the core where it is dissipated by the neutrinos. The 
heat diffusion through the crust can be described by the 
relativistic heat equation [111 ] : 
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where T is the temperature, t is the time, K is the ther- 
mal conductivity, CV is the specific heat and Q v is the 
neutrino emissivity. The effect of the gravity is given 
through the gravitational potential (f>, which enters in 
the definition of the readshifted temperature T = Te^, 

and the quantity T(r) = (l — 2Gm(r)/rc 2 ) , where 
G is the gravitational constant and m(r) is the gravita- 
tional mass included in a sphere of radius r. The latter is 
obtained solving the Tolman-Oppcnhcimcr-Volkoff equa- 
tion. 

In principle, the radial mass distribution and the prop- 
erties of the inner crust employed in the heat equation 
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should be determined consistently from the same equa- 
tion of state (EOS) of nuclear matter. To describe in an 
unified way both the crust and the core here we use the 
EOS given by the Skyrme interaction SLy4 [l§], which 
provides a reasonable star model [l3j] for the present 
study. The cooling calculations are done for a neutron 
star with a mass equal to 1.44 M Q . For this mass it is 
predicted a total star radius of 11.59 km and a central 
density of 3.51 po, where po — 2.9 x 10 14 g.cm -3 is the 
saturation density for symmetric nuclear matter. The in- 
ner crust, defined here as the part of the star with the 
density ranging between po/2 and pdrip — 3.285 x 10 11 
g.cm -3 , extends from i? c =10.74 km, which is the radius 
at the core-crust interface, to 11.29 km. 

The heat equation (1) is solved in the region of the 
inner crust matter. Since during the fast thermalisation 
the energy loss due to the neutrino emission in the crust 
is negligible we take Q u =0. The initial temperature dis- 
tribution in the crust and the boundary condition at the 
core-crust interface are chosen to simulate a rapid cooling 
process. Thus, at t = it is supposed that the crust has a 
flat temperature, i.e., T(r, t = 0) = 1}. We also consider 
that at the outer border of the crust the gradient of the 
temperature vanishes. 

The temperature evolution at the core-crust interface 
is obtained from Eq.(l) supposing that the radial distri- 
bution of the core temperature is uniform, dT/dr = 0. 
From Eq. (1) one thus gets 

(2) 



dT(r = R c ,t) 



dt C v 
Considering that the neutrino emissivity in the core is 
given by Q v = Q/T 9 6 [3, where T 9 = T/10 9 K, and sup- 
posing that during the crust thermalisation the specific 
heat remains constant in time we finally obtain 



C v 



T(r = R c , t) = % (l + eTft) 1/5 with e = 

(3) 

For the constant Qf, characterizing the fast neutrino 
emission in the core matter, we have taken the value 
Q f = 10 26 erg.cm-^s" 1 

The physical quantities which are essential for the crust 
thermalisation are the thermal conductivity and the spe- 
cific heat. In the inner crust the thermal conductivity 
is primarily determined by the electrons. Here we use 
the thermal conductivity parametrized by Lattimer et 
al[l| starting from the calculations of Itoh et al [HI]. For 
temperatures above 10 8 K, as used in this paper, the con- 
ductivity is nearly independent of the temperature and 
is given by K — C(p/ po) 2 ^ 3 , where C — 10 21 ergs cm -1 
s -1 and p is the baryonic density. 



III. SPECIFIC HEAT OF THE INNER CRUST 
MATTER 

The specific heat of the inner crust has contributions 
from the electrons, the lattice and the unbound neutrons. 



They are calculated for a given set of densities in the 
Wigner-Seitz approach. Since the electrons are ultrarel- 
ativistic, they are considered as an uniform degenerate 
gas with the specific heat given by [l6j] : 
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where V is the volume of the Wigner-Seitz cell and Z is 
the number of the electrons in the cell ( which is equal 
with the number of protons). 

For the specific heat of the lattice we use the approx- 
imation employed in Ref.0, i.e., C$ ttlce = 3k B /V, 
where, as above, V is the volume of the cell and ks is 
the Boltzman constant. 

We shall now discuss the temperature dependence of 
the specific heat of neutrons, which requires more elabo- 
rate calculations. The specific heat is calculated from the 
quasiparticle energies obtained solving the HFB equa- 
tions at finite-temperature in the Wigner-Seitz approxi- 
mation. The details of the HFB calculations in a Wigner- 
Seitz cell are given in Ref. [8j. In the HFB calculations 
the mean field is described with the same interaction used 
in the star model, i.e., the Skyrme force SLy4. A com- 
pletely consitent cooling simulation would require a inner 
crust with a Wigner-Seitz cells structure calculated from 
the same interaction. Since at present such cells structure 
is not available, here we shall use the Wigner-Seitz cells 
determined in Ref. jT8j . Their properties are summarize 
in Table I (see the Appendix). 

The pairing correlations in the inner crust matter are 
described with a density dependent contact force of the 
following form [l9j : 



V(r - r') = V [l - v (^-) a ]5(r - r'), 
Po 



(5) 



where p(r) is the baryonic density. To analyse the de- 
pendence of the crust thermalisation on the intensity of 
pairing correlations, in the calculations we have used two 
sets of parameters for the pairing force. They are cho- 
sen to simulate two possible scenarios for pairing in neu- 
tron matter corresponding to : (1) BCS calculations with 
realistic two-body interactions extracted from nucleon- 
nucleon scattering [||; in this approximation the max- 
imum pairing gap in uniform neutron matter is about 
3 MeV; (2) calculations which go beyond the BCS ap- 
proximation by taking into account in-medium effects 
on two-body interaction and self-energy; here we shall 
consider those calculations which predicts a maximum 
pairing gap in neutron matter of about 1 MeV [2(| Hl[ • 
These two scenarios, called below strong and weak pair- 
ing, can be simulated by using two pairing forces with the 
same parameters for the density dependent part, namely 
rj = 0.7 and a = 0.45, and two different strengths, i.c, 
Vb = {—570,-430} Mev /to -3 . These parameters have 
been used with an energy cut-off in the quasiparticle 
spectrum, required by the zero range of the pairing force. 
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FIG. 1: Neutron-specific heats in various Wigner-Seitz cells FIG. 2: Neutron-specific heats in various Wigner-Seitz cells 



for strong pairing. The specific heat is given in units of Boltz- 
man constant ks- 



for weak pairing. The specific heat is given in units of Boltz- 
man constant ks- 



The cut-off was introduced smoothly, i.e., by an exponen- 
tial factor e ~ B ' 2/10 ° acting for E % > 20 MeV, where E, 
are the HFB quasiparticle energies. 

With the setting discussed above we have solved the 
HFB equations for a given Wigner-Seitz cell and deter- 
mined the quasiparticle spectrum Ei and the correspond- 
ing entropy, i.e., 

S = -kg + !)(/< + (1 - fi) Ml - /*))• (6) 



where fi = [1 + exp(Ei / '/c^T)] -1 is the Fermi distribution 
and T is the temperature. From the entropy we then 
calculated the specific heat of the neutrons 



_TdS 



(7) 



where V is the volume of the Wigner-Seitz cell. 

The specific heat was calculated for all temperatures 
needed in the thermalisation process. Here we show the 



results for temperatures up to 500 keV. How looks the 
temperature dependence of specific heat for the two sce- 
narios of the pairing intensity is shown in Figs. (1,2). As 
noticed already in Ref.Q, we can see that the specific 
heat has very different values for the two pairing scenar- 
ios. It is also interesting to notice that the specific heats 
of the cells have rather different temperature dependence. 
Thus, for the strong pairing scenario, shown in Fig. 1, 
the specific heat is in the superfluid regime for the first 5 
cells (upper panel). This is not the case for the next cells 
(bottom pannel) where in the same temperature range 
there is a transition from the superfluid to the normal 
phase, as clearly seen for the cells nr. 6-8 (for the last 
two cells the transition temperature cannot be noticed 
because is too small). On the other hand, as seen in Fig. 
2, for the weak pairing the specific heat is entirely in the 
superfluid regime only for the first two cells. 

To illustrate the particular behavior of the specific heat 
in non uniform matter and the validity of various approx- 
imations, in what follows we shall discuss in more detail 
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FIG. 3: Neutron-specific heat in the the Wigner-Seitz cell 
nr 6 for strong pairing. The results corresponds to various 
approximations discussed in the paper. The specific heat is 
given in units of Boltzman constant ks ■ 



the results for the cell nr. 6, which contains N=460 neu- 
trons and Z=40 protons (see Appendix A). In this cell 
the HFB calculations predict 378 unbound neutrons. It 
is interesting that in spite of many neutrons in the cell 
the number of the bound neutrons in the cluster with 
Z=40 protons is equal to the magic number 82, as for the 
dripline nucleus 122 Zr (see, e.g., Ref.[22]). The specific 
heat given by the HFB spectrum, in which the contribu- 
tion of the cluster is included, is shown in Fig. 3 by full 
line. In the same figure are shown also the specific heats 
corresponding to two approximations employed in some 
studies P, 0] . In these approximations the non uniform 
distribution of the neutrons is replaced with a uniform 
gas formed by the total number of neutrons in the cell or 
by taking only the number of the unbound neutrons. The 
latter case is considered as an effective way of taking into 
account the influence of the cluster Q. How these ap- 
proximations work is seen in Fig. 3. To make the compar- 
ison meaningful, the calculations for the uniform neutron 
gas are done solving the HFB equations with the same 
boundary conditions as for the non uniform system, i.e., 
neutrons+cluster. As seen in Fig. 3, the transition from 
the superfluid to the normal phase is taking place at a 
lower temperature in the case of uniform neutron gas, es- 
pecially when are considered only the unbound neutrons. 
We can also notice that, in contrast to the uniform sys- 
tem, in the non uniform system the transition from the 
superfluid to the normal phase is smooth. 

To see better what happens in a non uniform system, 
in Fig. 4 is shown the evolution with the temperature of 
the pairing field in the cell nr. 6. We can notice that at 
zero temperature the pairing field is much larger in the 
surface of the cluster than in the bulk region. Due to this 
fact by increasing the temperature the pairs correspond- 
ing to the neutrons localized preferentially in the surface 



OF 1 I 1 I . _ . . !-j.-.-ali = =.-.'—■ ' — N — 




10 20 30 40 

R [fm] 



FIG. 4: Radial distribution of the pairing field for the Wigner- 
Seitz cell nr. 6 for various temperatures. The results corre- 
spond to a strong pairing. 
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FIG. 5: Neutron-specific heat in the Wigner-Seitz cell nr. 6 
for strong pairing. Here the results are shown up to high 
temperatures in order to illustrate the second transition re- 
gion around T= 900 keV. The specific heat is given in units 
of fcs. 



region of the cluster are destroyed gradually and much 
slower compared to the pairs formed by neutrons local- 
ized far from the surface of the cluster. In fact, as shown 
in Fig. 5, the non uniform system shows two transition re- 
gions, one around T=200 keV, corresponding mainly to 
the neutrons located far from the surface of the cluster, 
and another one, much less pronounced, around T=900 
keV, which corresponds to the neutrons localised in the 
surface region of the cluster. 
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FIG. 6: Time evolution of the readshifted temperature Te^ 
inside the inner crust for Ti=500 keV. The variable r repre- 
sents the radius of the star in km. 



IV. CRUST THERMALISATION 

In this section we discuss the crust thermalisation ob- 
tained by solving the heat equation (1) with the specific 
heats presented in the previous section. The time evolu- 
tion of the temperature in the inner crust is illustrated in 
Fig. 6 which shows the results obtained for an initial tem- 
perature Tj=500 keV. One can see that the crust ther- 
malisation is strongly enhanced by the presence of pairing 
correlations. 

Of physical interest is the time evolution of the so- 
called apparent surface temperature seen by a distant 
observer, Ti n f. For a non-magnetized envelope composed 
by iron an d lig ht element the surface temperature T in f 
is given by [231 ] 



eff 



2GM 



(8) 



where 



10 b 



10 



14 



T eff 

C = T b9 - (T*/10 



( (7C) 2.25 + (c/3) l 



(1/4) 



g = GM/R^l - GM/c 2 R 
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In the equations above T» = (7T b9y /gi4) " where gn 
g/10 14 , T b9 = T b /10 9 and T b is the temperature at the 
outer edge of the inner crust provided by the heat equa- 
tion. 

The time evolution of the surface temperature Ti„f is 
displayed in Fig. 7 for three initial temperatures of the 
crust, Ti = {100,300,500} keV. We can clearly notice 
that the pairing is reducing significantly the thermalisa- 
tion time of the crust. For example, as seen in Fig. 7, bot- 
tom panel, the cooling time drops from about 75 years, 
the value obtained when the neutrons are considered in 
the normal phase, to about 55 years for weak pairing 
and to about 30 years for strong pairing. Thus the two 
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FIG. 7: Time evolution of the surface temperature for the 
initial temperatures Ti = {100, 300, 500} keV. 



pairing scenarios give very different predictions for the 
thermalisation time. The same conclusion was obtained 
previously in a schematic cooling model Q • However, the 
latter predicts very different thermalisation times com- 
pared to the realistic cooling model employed here. 

In Fig. 7 are shown also the thermalisation times ob- 
tained by neglecting the effect of the clusters, i.e., sup- 
posing that the baryonic matter in the inner crust is uni- 
form. In this case the specific heat of the neutrons is 
calculated from the quasiparticle spectrum of BCS equa- 
tions solved for infinite matter. For the BCS calculations 
we have used the single-particle energies provided by the 
asymmetric nuclear matter with the same neutron to pro- 
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ton fraction as in the Wigne-Seitz cells of Table I. As seen 
in Fig. 7, the results for the weak pairing scenario indi- 
cate that the surface temperature is dropping faster for 
uniform matter than for non uniform matter. For the 
strong pairing the behaviour is opposite. Thus the clus- 
ters have a non-negligible and a non-trivial effect on the 
time evolution of the surface temperature. 

In order to test the effect of the protons on uniform 
matter calculations we have also evaluated the time evo- 
lution of the surface temperature considering only a uni- 
form gas of neutrons, i.e., neglecting in the BCS equa- 
tions the contribution of the protons to the single-particle 
energies of neutrons. The results are very similar to the 
ones obtained for the uniform asymmetric matter. 

V. SUMMARY AND CONCLUSIONS 

In this paper we have studied how the thermalisation of 
neutron stars crust depends on pairing properties and the 
cluster structure of the inner crust matter. The thermal 
evolution was obtained by solving the relativistic heat 
equation with initial conditions specific to a rapid cool- 
ing process and supposing that during the thermalisation 
there are no sinks or sources of energy in the crust. The 
specific heat of neutrons was calculated from the HFB 
spectrum, taking into accound consistently the effects of 
nuclear clusters, pairing correlations and temperature. 
The thermal evolution of the inner crust was analysed 
using for the neutrons two sets of specific heats obtained 
with a strong and a weak pairing force which simulate 
two possible scenarios for the intensity of pairing correla- 
tions in neutron matter. The results show that the crust 
thermalisation is strongly influenced by the intensity of 
pairing correlation. This result confirm what it was found 
earlier in Ref. [1] with a schematic cooling model. How- 
ever, the latter predicts thermalisation times which are 
very different from the results obtained with the realis- 
tic cooling model employed in this study. We have also 



shown that the cluster structure of the inner crust affects 
significantly and in a non-trivial way the time evolution 
of the surface temperature, mainly for weak pairing and 
before the thermal equilibrium between the crust and the 
core is reached. 
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VI. APPENDIX A 

In this appendix are summarized the properties of the 
Wigner-Seitz cells determined in Ref. [18| and used in 
this paper. Compared to Ref. [hsl]. here we didn't include 
the most dense cell, which is close to the region of more 
complicated pasta phases. 

cell N Z Rws p 
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[g.cm- 3 ] 
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140 


40 
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4.7 x 10 11 
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160 


40 


49 


6.7 x 10 11 
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210 


40 


46 


1.0 x 10 12 
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280 


40 


44 


1.5 x 10 12 
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460 


40 


42 


2.7 x 10 12 


5 


900 


50 


39 


6.2 x 10 12 


4 


1050 


50 


36 


9.7 x 10 12 


3 


1300 


50 


33 


1.5 x 10 13 


2 


1750 


50 


28 


3.3 x 10 13 


1 


1460 


40 


20 


7.8 x 10 13 



TABLE I: The structure of the cells determind in Ref. [3 , 
i.e., the baryonic densities (p), the number of neutrons (N), 
the number of protons (Z) and the cell radii (Rws) 
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